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Abstract 

We ascertain the effectiveness of the second post-Newtonian approximation 
to the gravitational waves emitted during the adiabatic inspiral of a compact 
binary system as templates for signal searches with kilometer-scale interfer- 
ometric detectors. The reference signal is obtained by solving the Teukolsky 
equation for a small mass moving on a circular orbit around a large nonrotat- 
ing black hole. Fitting factors computed from this signal and these templates, 
for various types of binary systems, are all above the 90% mark. According 
to Apostolatos' criterion, second post-Newtonian waveforms should make ac- 
ceptably effective search templates. 

Pacs numbers: 04.25.Nx, 04.30.Db, 04.80.Nn 
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I. INTRODUCTION 



The current construction of kilometer-scale interferometric gravitational- wave detectors, 
such as the American LIGO (Laser Interferometer Gravitational- wave Observatory U) and 
the French- Italian VIRGO @, has motivated a lot of recent work on data analysis tools 
to search and measure the expected gravitational- wave signals 0. Thus far, the mostly 
studied waves have been those generated by the adiabatic inspiral of a compact binary 
system, composed of neutron stars and/or black holes ||. Because of their reliably large 
event rate ^-0], these waves are currently believed to be the most promising for detection 
by LIGO/VIRGO. And because their waveforms can be accurately predicted, searching and 
measuring these waves will be best accomplished with the well-known technique of matched 
filtering ||. In this paper we focus on certain issues related to the search of gravitational 
waves from inspiraling compact binaries. 

The method of matched filtering || relies on the fact that an accurate model can be 
constructed for the expected signal. The model (usually called the template) depends on 
various parameters, such as the signal's time of arrival, its initial phase, the masses of the 
companions, their spin angular momenta, and others. The template is cross-correlated with 
the detector output, and the parameters are varied until the largest correlation is obtained. 
If the maximum correlation (usually called the signal-to-noise ratio) is larger than a specified 
threshold value, then a signal is concluded to be present. 

A large effort is currently underway to calculate inspiraling binary waveforms to high 
order in a post-Newtonian approximation p,[10| . This approach is based on the assumption 
that the orbital motion is sufficiently slow that the waves can be expressed as a power series 
in v/c, the ratio of the orbital velocity to the speed of light. Because the series converges 
very slowly, if at all, a large number of terms are required [[Ll|,[rj|. Thus far, the waveforms 
have been computed to order (v/c) 5 (2.5 post-Newtonian, or 2.5PN, order) beyond the 
leading-order, quadrupole-formula expressions [T(| . 

Central to the theory of matched filtering is Wiener's theorem ||, which states that 
the expectation value of the correlation is maximized when the template parameters match 
those of the actual signal. For the theorem to apply, however, the functional form of the 
template must be the same as that of the signal, so that template and signal are allowed 
to disagree only in the value of their parameters. In view of the fact that the signal is 
governed by the exact laws of general relativity, while the template necessarily constitutes 
an approximation, we must expect that the actual signal-to-noise ratio will be somewhat 



less than the maximum allowed by Wiener's theorem [13]. We must also expect that the 



value of the template parameters which maximize the correlation will not precisely match 
the value of the source parameters. 

Apostolatos |14] has introduced the fitting factor (FF) as a measure of template imper- 
fection. The fitting factor is the ratio of the actual signal-to-noise ratio, obtained with the 
imperfect template, to the signal-to-noise ratio that would be obtained if a perfect template 
were available [|T^]. A fitting factor of unity means that the template is a very accurate 
representation of the signal (a perfect template). A value less than unity means that the 
template reproduces the signal only imperfectly. The loss in event rate due to template 
imperfection is given by 1 - FF 3 0. 

The fitting factor is computed as follows jl4jl . Let s(t) be the actual signal, and let h{t; 6) 
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be the templates, with the vector 6 representing the parameters. We denote the Fourier 
transforms of these functions by s(f) and h(f;0), respectively, where a(f) = J a{t)e 2m ^ dt 
for any function a(t). We define the ambiguity function A(0) by 



m = (i) 

\s\s){h\h) 



where, for any functions a(t), b(t), 



oo 



m = 2 , HfED_tWl d , (2) 

Jo S n (f) 

Here, an overbar indicates complex conjugation, and S n (f) is the (one-sided) spectral den- 
sity of the detector noise. For interferometric detectors of the advanced-LIGO type, it is 



appropriate @ to set S n (f) = S [(/ min //) 4 +2+2(/// min ) 2 ] for / > 10 Hz, while S n (f) = oo 
for / > 10 Hz; Sq is a normalization constant irrelevant for our purposes, and f m \ n = 70 Hz 
is the frequency at which the noise is minimum. The quantity {a\b) is usually interpreted as 
an inner product in the Hilbert space of gravitational-wave signals |T5[ , and arccos^4(0) is 
the angle between the signal s(t) and the template h(t; 0). The fitting factor is given by 

FF = max„4(0), (3) 
and arccosFF is the minimum value of this angle. Apostolatos has suggested that 



templates for which FF > 90% would be acceptably accurate, leading to a loss in event rate 
no larger than approximately 27%. 

To compute the fitting factor one must first provide a reference signal s(t) and a param- 
eterized set of templates h(t; 0). While post-Newtonian waveforms make an obvious choice 
of templates, the selection of a reference signal is much more delicate, because an exact 
representation of the gravitational waves is not available. In all studies undertaken thus 
far [|14] , |T6| - 21], the reference signal was also chosen to be a post-Newtonian waveform, but 



one computed with more accuracy than the templates. (For example, the signal is taken to 
be a 2.5PN waveform, while the templates are chosen to be Newtonian, or 1PN, or 1.5PN, 
or 2PN waveforms.) This approach is flawed, because even a high-order post-Newtonian 



approximation may still be a poor representation of the true general-relativistic signal |]TT 

In this paper we make a different choice of reference signal which, we believe, bears a 
closer resemblance to the true general-relativistic signal for binary systems with small mass 
ratios. This signal is obtained by solving the Teukolsky equation |22j for the gravitational 
perturbations produced by the circular motion of a small mass around a large nonrotating 
black hole fl2"3|j2"4jl . In this approach, no assumption is made regarding the size of v, the orbital 
velocity. This is one advantage over post-Newtonian theory. However, one disadvantage is 
the necessary restriction to binary systems with small mass ratios. 

Our fitting factors are therefore computed by taking s(t) to be derived from the equations 
of black-hole perturbation theory, and by taking h(t; 6) to be given by the so-called restricted 
second post-Newtonian approximation. (This terminology will be explained below; we use 
2PN waveforms because the formally more accurate 2.5PN waveforms are known to be a 
poorer representation of the true signal fll|| .) We consider the binary systems listed in 
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TABLE I. The binary systems considered in this paper are listed in the increasing order of 
their total mass. The individual masses mi and 7712, and the chirp mass A4 a ctuai are given in solar 
units. The mass-ratio parameter r/ actua i is dimensionless. All quantities are defined in the text. 

System mi m 2 A4 actU ai factual 



1.4 + 1.4 
0.5 + 5 
1.4 + 10 
10 + 10 
4 + 30 



1.4 
0.5 
1.4 
10.0 
4.0 



1.4 
5.0 
10.0 
10.0 
30.0 



1.2188 
1.2322 
2.9943 
8.7055 
8.7340 



0.2500 
0.0826 
0.1077 
0.2500 
0.1038 



TABLE II. Searching perturbation-theory signals with second post-Newtonian templates. For 
each binary system, and for two types of signals, we list the values of A4* , rj*, and t* corresponding 
to the global maximum of A(t* , A4* , rf). The fitting factor is also listed. All quantities are defined 
in the text. For all numbers the numerical uncertainty is ±5 in the last digit. 

System Signal M* max r/^ ax ^ ax FF 



1.4 + 1.4 


one-mode 


1.0033 


0.401 


0.3608 x 10 -4 


94.7% 


1.4 + 1.4 


I < 3 


1.0033 


0.402 


0.3562 x 10~ 4 


93.0% 


0.5 + 5 


one-mode 


1.0017 


0.550 


0.2345 x 10~ 4 


97.8% 


0.5 + 5 


I < 3 


1.0017 


0.550 


0.2421 x 10" 4 


95.2% 


1.4 + 10 


one-mode 


1.0013 


0.496 


-0.2186 x 10~ 3 


95.5% 


1.4 + 10 


I < 3 


1.0013 


0.496 


-0.2154 x 10~ 3 


91.4% 


10 + 10 


one-mode 


1.0001 


0.388 


-0.1197 x 10~ 3 


98.4% 


10 + 10 


I < 3 


1.0001 


0.388 


-0.1179 x 10" 2 


91.7% 


4 + 30 


one-mode 


0.9859 


0.367 


-0.5256 x 10~ 2 


95.4% 


4 + 30 


I < 3 


0.9859 


0.367 


-0.5192 x 10~ 2 


86.6% 



Table I, which consist of three low mass-ratio systems (O.5M + 5M , 1.4M + 10M Q , and 
4M Q + 3OM , where M is the solar mass) and two equal-mass systems (1.4M + 1.4M 
and 1OM + 1OM ). Although our analysis is not expected to be valid for the equal- mass 
systems, we nevertheless include them as an indication of the robustness of our conclusions. 
For all binary systems the companions are assumed to be nonrotating. 

Our results, displayed in Table II, show that the 2PN waveforms make acceptably accu- 
rate templates, in the sense of Apostolatos' criterion [14 1 : except for the large-mass system 
4M + 3OM , the fitting factors are all above the 90% mark. This is the main conclusion of 
this paper. We must point out that our analysis ignores such important issues as discrete 
template spacing p5|j26|] , computing power, and spin- induced signal modulations fi~4| , |2~7 . 
We leave such considerations for future work. Throughout the paper we work in units such 
that G = c = 1. 
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II. REFERENCE SIGNAL 



The formalism of black-hole perturbation theory [E3J returns the following expression for 



the gravitational waves emitted by a mass m\ in a fixed circular orbit around a nonrotating 
black hole of mass m 2 ^> m\. 



S + (t) - lS x (t) (X^2 A l™{v)e~ 



(4) 



1=2 m=-l 



Here, s + (t) and s x {t) are the two fundamental polarizations of the gravitational waves, 
Ai m {v) are mode amplitudes, and -2Yi m (9,(fi) are spin-weighted spherical harmonics p8| . 
with 9 and the polar angles of the gravitational-wave detector. The missing normalization 
constant is unimportant for our purposes. In Eq. (g), the orbital velocity v and the orbital 
angular velocity Q are related by v 3 = MQ, where M = m\ + m 2 is the total mass. The 
mode amplitudes Ai m (v) must be calculated separately for each value of v by numerically 
integrating the Teukolsky equation 

Equation (g) describes the waves produced by a mass moving on a fixed orbit. This 
expression is easily generalized to describe an adiabatic inspiral: the orbital velocity v now 
becomes a function of time, which is determined by 



dv 
~dt 



32^ g P(v) 



5M 2 Q(v) 



(5) 



Here, /i = mim 2 /M is the reduced mass, P(v) = (dE / dt) / (dE / dt) N is the rate at which the 
gravitational waves remove orbital energy from the system, normalized to the quadrupole- 
formula expression (dE/dt) N = -(32/5)(>/M) V°, and Q{y) = (dE/dv)/(dE/dv) N = (1 - 
6V 2 ) (1 — 3t> 2 )~ 3 / 2 gives the differential relation between orbital energy and orbital velocity, 
normalized to the Newtonian expression (dE/dv)N = —y/v. The function P(v) must be 
computed numerically |24 ]. As a consequence of the time dependence of Q = v 3 /M, the 
factor nt in Eq. (g) must be replaced by $(t) = /* fi(t') dt'. 

To simplify the expressions we assume that the detector is situated at (p = 0, so that the 
spherical harmonics are functions of 9 only. (Because the orbit is circular, this assumption 
implies no loss of generality.) We also assume that the detector is oriented in such a way that 
it is sensitive only to the "+" polarization Simple manipulations, using the relation 

Ai_ m = (-l) l Ai m 0, then bring Eq. (g) to the form 



s(t) oc £ £ Re 

1=2 m=l 



where s(t) = s + (t), 



A lm (v)e- im *W] S lm (6) 



(6) 



fmdS lm (6) = -. 2 Y lm (e, 0)- 
by Eq. (|). 



. 5M rv Q{v') , , 
V ; 32/1 J v' 6 P(v') ' V ; 

'—l) l _2Yi- m {9, 0). The relation between t and t> is given implicitly 



5 



Equations (|6]) and (^) give the reference signal in the time domain. Computation of the 
fitting factor, however, requires an expression in the frequency domain. We must therefore 
Fourier transform Eq. (|6]), which we do in the stationary phase approximation. Skipping all 
intermediary steps, we find 



where 



1=2 m=l 



Q(v) 



nl/2 



mv ll P(v) 



A lm (v)e im ^S lm (9), 



(8) 



5M 
32^ 



(v d - v'^Qjv') 
v' 9 P(v') 



dv . 



(9) 



In these expressions, / and v are related by v — (27rM//m) 1//3 , and vq is an arbitrary 
constant which sets the zero of the phase function \l/(f), which must also be evaluated 
numerically. 

In Eq. (pi), the sum over / must necessarily be truncated at some finite value. The 
dominant part of the signal comes from the I = m = 2 mode. We shall refer to a signal 
containing only this mode as a one-mode signal. The second largest contribution to the 
signal comes from the I = m = 3 mode, and in fact, all additional contributions can be 
neglected. We shall therefore truncate the sum at I — 3, having numerically verified that 
additional terms are indeed irrelevant. Because S33/S22 = — (21/10) 1 / 2 sin 9, the relative 
importance of the I = 3 terms is maximized when 9 = ir/2. We shall therefore set 9 to this 
value throughout the paper. 

The signal must be cut off when the orbiting mass reaches the innermost stable circular 
orbit at v = Because the I = m = 2 mode dominates the signal, we simply let 

§(/) = when / > (6 3 / 2 7rM) -1 . 

III. TEMPLATES 



The restricted second post-Newtonian approximation to the gravitational-wave signal is 
given in the frequency domain by 



h(f;e)<xf- T ' b e 



7/6 JHf) 



(10) 



where = 2vr/t c - C + <f>(f) 

3 



128 



(nMf) 



-5/3 



20 /743 11 

1 H 1 v ) x" 

9 V336 4 ' ' 



10 



3058673 5429 



617 



1016064 



-7? H 77" 

1008 ' 144 ' 



x 



(11) 



with x 3 = irr]~ 3 / 5 Aif '. Here, {4> c , t c , Ai, rj} are the template parameters, with <p c the phase 
at coalescence (formally / = 00), t c the time at coalescence, A4 the chirp mass, and r\ the 
mass-ratio parameter. The term "restricted" refers to the fact that while the phase function 
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<f>(f) is calculated to 2PN order beyond the quadrupole-formula expression, the amplitude 
function /™ 7//6 is kept at leading order. 

Because the templates do not have the same functional form as the signal, M. and r\ do 
not have direct physical meaning. However, if the templates were an accurate description 
of the waves, then the signal would be best reproduced when M. = M. ac t U ai and 77 = factual, 



where 31 



(m 1 m 2 ) 3/5 m 1 m 2 , , 

-^actual / \i /c j factual / \ 9 j (-^J 

(mi + m 2 ) i° {rrii + m 2 ) 

with mi and 7712 denoting the actual masses of the companions. It is therefore convenient 
to introduce the rescaled parameters A4* and rj*, defined by 

M = M*M acWal , 1] = r/*r/ actual . (13) 

It is also convenient to make the following transformations: 

c = o + 0*, t c = t (l + t*), (14) 

where 

<j>o = 2irf t c + <j>(f ), *o = -^^o)« ( 15 ) 

with fo an arbitrarily chosen frequency. Their purpose is to ensure that when cj>* = t* = 0, 
then V>(/o) — dip/df(f ) = 0. These relations are analogous to those holding for ty(v) at 
v = v [see Eq. (0)]. We shall henceforth take = (<f>*,t*,M*,7)*) to be our template 
parameters. 



IV. COMPUTING THE FITTING FACTOR 

Evaluation of the ambiguity function, Eq. fll]), for a given binary system and a given 
set of template parameters, is a straightforward numerical problem. First Ai m (v) and P(v) 
are computed once and for all, and tabulated for many values of the orbital velocity. Next, 
*i?(v) is calculated for the binary system under consideration. Combining the results and 
summing over the relevant modes gives s(f), Eq. (^). On the other hand, h(f; 6) is given 
in analytical form by Eqs. (|10[) — (|T5|) . The overlap integral (s\h) is then evaluated. To 
do this efficiently it is advantageous to choose the constants fo and vq [see Eqs. (P) and 



(15)] so that ^/(/), "4>(f), and their slopes all vanish where the integrand is largest. This 
ensures that signal and template are approximately in phase in the frequency interval that 
contributes the most to the integral. Now, the integrand is fairly well approximated by 
\h\ 2 /S n (f) oc z~ 7 / 3 /(z~ 4 + 2 + 2z 2 ), where z = ///mm, which is maximum at z ~ 0.6654. 
We therefore set fo = 46.58 Hz. And since s(f) is dominated by the I = m = 2 mode, we 
correspondingly put t>o = (7rM/ ) 1//3 . With these choices, computation of (s\h) is optimized. 
Evaluation of (s\s) and (h\h) presents no additional difficulties, and combining the results 
gives A{0). 
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FIG. 1. Shown here are surfaces of constant ambiguity function A(t* , M* , 77*) for the binary 
system 1.4M + IOMq. The reference signal is a one- mode signal, and the templates are second 
post-Newtonian waveforms. The light grey surfaces represent A = 15%, the dark grey surfaces 
represent A = 30%, and the black surfaces represent A = 50%. The global maximum of 95.5% is 
contained within the large structure to the left of the box. The other structures contain additional 
local maxima, and these make an automated search for the true maximum difficult. The region of 
parameter space covered by the box is -0.06 < t* < 0.35, 0.96 < M* < 1.05, and 0.35 < rf < 0.90. 



We now discuss the issue of maximizing the ambiguity function over the template pa- 
rameters. Maximization over 0* is straightforward. Since 0* comes as an overall phase 
in h(f;0), we have that {s\h) = Re(|/|e ia ), where a is linear in 0*, and where I is some 
function of the remaining parameters. Maximization over 0* is achieved by making a = 0, 
which yields (s\h) = Since (h\h) and (s\s) do not depend on 0*, we have also maximized 
A{6) over 0*. We will denote this reduced ambiguity function by A(t*, Ai*, r]*). Maximizing 
this over the remaining parameters is not a straightforward numerical problem, because the 
ambiguity function displays a very rich structure, featuring many local minima and maxima 
(see Fig. 1). We proceed as follows. For each binary system under consideration we evaluate 
the reduced ambiguity function in a fairly large neighborhood of (£*, A4*, 77*) = (0, 1, 1). The 
function is then displayed graphically with the help of a three-dimensional visualization soft- 
ware [52 1, and the approximate position of the global maximum is determined. Finally, the 
maximum is determined accurately by applying Powell's method as described in the book 
Numerical Recipes [|33|]. The value of A(t* , A4* , i]*) at the global maximum is the fitting 
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TABLE III. Searching a one-mode signal with Newtonian templates. For each binary system 
we list the values of Ai* and t* corresponding to the global maximum of A(t* , Ai*). The fitting 



factor is also listed. The numerical uncertainties are as in Table II. 

System M* nax 



1.4 + 1.4 0.9944 -0.4120 x 10" 3 79.2% 

0.5 + 5 0.9999 -0.3597 x 10" 4 51.6% 

1.4 + 10 1.0201 -0.1057 x 10~ 2 55.8% 

10 + 10 1.0563 -0.1361 x lO" 1 70.1% 

4 + 30 1.1503 -0.2339 x 10- 1 61.3% 



factor FF. 

The calculations described in this section require a great deal of numerical work. To 
minimize the possibility of making coding errors which would affect our results and conclu- 
sions, a separate code was written by each of the authors. Because the outputs of the two 
separate codes agree within the numerical uncertainty quoted in the caption of Table II, we 
are quite confident in the accuracy of our results. 



V. RESULTS AND DISCUSSION 

Our results are displayed in Table II. For each binary system we consider two types of 
signals. The one-mode signal contains only the I — m — 2 mode, while the I < 3 signal 
contains all modes up to I = 3. Inclusion of additional modes does not change the results 
within the numerical uncertainty. We see that the fitting factors for the one-mode signals 
are large, all above approximately 95%. However, inclusion of additional modes decrease 
the fitting factors appreciably: the reduction ranges from approximately 3% for low-mass 
systems to approximately 7% for large-mass systems. We have verified that most of this 
reduction can be attributed to the I = m = 3 mode. Nevertheless, the fitting factors are all 
larger than 90%, except for the large-mass system 4M + 3OM . The second post-Newtonian 



templates therefore meet the Apostolatos criterion | 14| , and we conclude that they should 
be adequate to search for gravitational waves emitted during the adiabatic inspiral of a 
compact binary system. 

For comparison, we show in Table III the results obtained when Newtonian templates 
are used. [We discard all post-Newtonian corrections in Eq. (|i~l~D, and rf disappears from 
the list of template parameters.] Here the fitting factors are computed using a one-mode 
signal, as this gives the largest result. We see that these fitting factors are all much smaller 



than those of Table II. The Newtonian templates do not meet the Apostolatos criterion |14 

The reduction in fitting factor that occurs when a one-mode signal is replaced by a / < 3 
signal is clearly caused by the template's inability to keep track of the additional frequency 
components. We have attempted to improve the performance of our templates by adding a 
term 

^7/3 

h ^f. 0) oc ^_ (1 _ 4^)1/2 f-V* x sin0e i[^/tc-3«c/2+^2//3)] j (16) 
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where the constant of proportionality is the same as in Eq. (fL0|) . While h(f;0) describes 
waves at twice the orbital frequency and is the leading-order term in a post-Newtonian 
expansion, h 3 (f; 0) describes waves at three times the orbital frequency and is a correction 
term of order v. Equation (|l^) was obtained by Fourier transforming Eq. (3b) of Ref. ||. 
Our expectation was that this additional term would do a good job at reproducing the 
behavior of the I = m = 3 mode, and that the loss in FF would be regained. Calculating 
fitting factors with h + h 3 as templates shows otherwise. (Some thought must be given to 
the fact that <fi* no longer appears as an overall phase in the templates.) The improvement 
in FF was less than 1% for the large-mass systems, and less than 0.1% for the low- mass 
systems. This prompts us to conclude that it would be pointless to introduce additional 
frequency components in the search templates. 

It is interesting to observe that for low-mass systems, the maximum of the ambigu- 
ity function is obtained when M. is approximately 0.1% above its actual value, while 77 is 
approximately 40% of its actual value. These numbers can be loosely interpreted as the sys- 
tematic errors incurred when attempting to estimate the source parameters with imperfect 
2PN templates. The systematic errors should be compared with the anticipated statistical 
uncertainties associated with the measurement of a noisy signal [|l^,^l|,^,0] . Cutler and 



Flanagan [TJ|] have determined that for low-mass systems known to be composed of nonro- 
tating companions, the statistical uncertainty associated with the estimation of M. is of the 
order of 0.005%, while it is of the order of 1% for r\ [0. (This assumes detection with a 
signal-to-noise ratio of 10.) Second post-Newtonian templates therefore give rise to system- 
atic errors that are much larger than the statistical uncertainties. While these templates 
could be adequate to search for signals, they should not be used as measurement templates. 
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